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ABSTRACT 



Stars within 0.1 pc of the supermassive black hole Sgr A* at the Galactic centre are 
expected to encounter a cluster of stellar-mass black holes (BHs) that have segregated 
to that region. Some of these stars will scatter off an orbiting BH and be kicked 
out of the Galactic centre with velocities up to ~ 2000 km s^^. We calculate the 
resulting ejection rate of hypervelocity stars (HVSs) by this process under a variety of 
assumptions, and find it to be comparab le to the tidal disruption rate of binary stars 
by Sgr A*, first discussed bv lHilld (|1988D . Under some conditions, this novel process is 
sufficient to account for all of the hypervelocity B-stars observed in the halo, and may 
dominate the production rate of all HVSs with lifetimes much less than the relaxation 
time-scale at a distance ~ 2 pc from Sgr A* (> 2 Gyr). Since HVSs are produced by at 
least two unavoidable processes, the statistics of HVSs could reveal bimodal velocity 
and mass distributions, and can constrain the distribution of BHs and stars in the 
innermost 0.1 pc around Sgr A*. 
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1 INTRODUCTION 

Hypervelocity stars (HVSs) have velocities so great (~ 
1000 kms^'^) that they are gravitationally unbound to the 
Milky Way galaxy. Since the discovery of the first HVS by 
iBrown et aiTll2005l). six additional HVSs have been found 
tedelma nn et all boosl : iHirsch et all l2005l : iBrown et all 
[2006a, b). The age and radial velocities of all but one HVS 
are consistent with them originating from the Galactic cen- 
tr43) the most natural site for producing them (Hills 1988). 
As more HVSs are found, they can be used to constrain 
many properties of both the Galactic centre as well as the 
Milky- Way galaxy as a whole, providing information about 
the G alactic potential (jGnedin et al.l l2005l : iBromlev et al.l 
20061). the merger his tory of the Sgr A* (jLevinl l2005l : 



the Galactic centre. iHiUsI I 1988h showed that when a tight 
stellar binary (with a separation < 0.1 AU) gets sufficiently 
close to a SMBH, it will be disrupted by its strong tidal 
field. Consequently, one member of the binary could be 
ejected from the Galactic centre with sufficient ener gy to 
escape the gravitational pote ntial of the entire galaxy (|Hillsl 
I1991I : iGualandris et al.|[2005l ). The other binary member is 
expected to rema in in a highly eccentri c orbit around the 



SMBH (see, e.g., iGould fc Quillenll20o3 : iGinsburg fc Loebl 



Baumgardt et al. I2OO6I ). as well as the mass distribution of 



stars near Sgr A*. 

Nearly 17 ye ars be fore the initial discovery by 

iBrown et all l|2005ll . iHillsl (|l988l ) first proposed that HVSs 
should populate the galaxy and provide indirect evidence 
for the existence of a supermassive black hole (SMBH) in 
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In the case of Edclma nn et al ] l|2005l ). the star appears to orig- 
inate from the Large Magellanic Cloud. 



'200d). lYu fc Tremainel (1200311 I analysed the production rate 
of HVSs in more detail, specifically for stars originating near 
the SMBH in our galaxy, Sgr A*. The a uthors also cor- 
rected the initial calculation of iHillsl { 19881 ) and accounted 
for the diffusion of hard binaries into the "loss-cone", find- 
ing the production rate of HVSs to be ~ 10~^ yr"^, nearly 
three orders of magn itude below the previous estimate. 
IYu fc Tremain^ l|2003l ) also looked at two additional mech- 
anisms for producing HVSs. They found that the rate could 
easily be higher (~ 10~^y r~^) if Sgr A* had a massive 
binary companion (see also iLevinI l2005l : iBaumgardt et aP 
l2006l : ISesana et al.ll2006l ) , and also determined that star-star 
scattering near Sgr A* resulted in a nearly undetectable rate, 
due to physical collisions among the two stars. 

In examining the rate estimate from IYu fc Tremainel 
(l2003l ) one may adopt a simple mass function (MF) for 
stars, An/Am oc to estimate the total number of B- 
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type HVSs in the Galactic halo. Using the Salpeter MF 
slope /3 = 2.35 and a minimum mass limit of 0.5 Mq one 
expects there to be about ~ 80 HVSs with masses be- 
tween 3 a nd 5 , consistent with the observed rate for 
such stars (|Brown et al.|[2006al ). However this is an overes- 
timate. The Salpeter MF overestimates the total number 
of stars with mass > 3Mq, since star formation occurs 
continuously over time near the Galactic centre and mas- 
sive stars have short lifetimes (| Alexander fc Sternberll 19991 : 
IPerets et al.ll2006l ). Accounting for continuous star forma- 
tion a s well as a more rea listic distribution of binary param- 
eters, iPerets et al] (|2006l ) found a lower total HVS rate of 
5 X 10~'^yr~^ yielding ~ 1 HVS in the entire galaxy with 
mass between 3 and 5 Mq . The discrepancy between this cal- 
culation and the number of observed HVSs may be overcome 
by massive perturbers in the Galactic centre, which reduce 
the relaxatio n time and increase the rate to that observed 
iPerets et aLll2006D . Nevertheless, their results are highly 
sensitive to the mass distribution of the perturbers. In their 
calculations the authors use a very flat slope for the distribu- 
tion of masses which has not been well constrained by obser- 
vations. In their estimate, they assume that all stars are ini- 
tially relaxed and therefore fill the entire energ y-momentum 
space. The lifetimes of the stars observed bv iBrown et al.l 
(|2006al ) are very similar to the enhan ced relaxation time , 
except in the most optimistic case of IPerets et all l|2006t) . 
W e argue then t hat th e total number of B-type HVSs given 
bv IPerets et al.l (|2006| ) may still be an overestimate, and an 
additional source of HVSs may be required. 

In this paper, we propose a novel source of HVSs in 
the Milky Way: in the dense stellar cusp near Sgr A* stars 
scatte r o ff stellar-m ass black holes (BHs) that are segregated 
there (|Morrislll993l ) and recoil out of the Galactic centre. The 
existence of stars on radial orbits o riginating from the Galac- 
tic ce ntre was first proposed by iMiralda-Escude fc Gouldl 
(|2000l ) as evidence for a segregated BH cluster. Here, we 
show that there should exist a high-velocity tail of stars, 
significant enough in nu mber to account for s ome if not all 
of the HVSs observed bv lBrown eTeI\ (12006^ ). 

This paper is organised as follows. In JHwe describe the 
distribution of stars and BHs around Sgr A*. We describe in 
detail our assumptions and calculations in 321 and present 
our results in Finally, in we discuss the observational 
implications and prospects for future work. 



2 THE CUSP OF STARS AND BHS AROUND 
Sgr A* 



Using near-infrared adaptive optics imaging, iGenzel et al.l 
l|2003t ) showed that the density of stars in the innermost 
~ 2 pc of the Galactic centre is well-fit by a broken power- 
law radial profile, p{R) = 1.2 x 10^(ii/0.4pc)-" M© pc"^ 
where a ~ 1.4 for R < 0.4 pc and a ~ 2.0 for R > 0.4 pc, as- 
sumi n g that Sgr A* is at a distance of 8 kpc (|Ghez et al.l 
I2OO5I : lEisenhauer et al.l |2005| ). Additional observations of 
the background light density show that the cusp might be 
slight ly shallower in the inner few arcseconds ( Schodel et al. 
I2OO7I ). Indeed, for a dynamically relaxed stellar system, 
stars near a SMBH are expected to be in a cusp pro- 
file with 1.5 < a < 1.75, depending on the ma ss distri- 
bution of the stars (|Bahcall fc Wolf Il97d 1 19771 : but see 



also lHopman fc AlexandedbOOebl l. lGenzel et~all (120031 ) nor- 
malised the density profile by assuming that the total mass 
of stars within 1.9 pc of Sgr A* is 3.2 x 10*^ Mq, and that 
the mass distribution follows the same density profile as the 
observed stars with a constant mass-to-light ratio. Invari- 
ably, stellar evolution and mass segre gation should alter the 
mass- to-light ratio at different radii l|Hopman fc Alexanderl 
l2006bl ). and cause the actual density profile to deviate from 
the observed number counts of stars. In addition to deter- 
mining that the stars are in a cusp profile, the observations 
also indicate that the stars ha ve formed mostly continuously 
with a standard i nitial MF ([Alexander fc Sternberg! Il999l : 
iGenzel et aLll2003l ). 

Besides the observed cusp of stars, there are two ad- 
ditional structures found near Sgr A* that suggest ongoing 
star formation very close to the SMBH. Within 1 arcsec 
(~ 0.04 pc) of Sgr A* is a nearly isotrop ic cluster of mas- 
sive B-type stars, the so -called "S-stars" (|Eisenhauer et ah! 
I2OO5I: iGhez et al.ll2005h . whose origin remain s a mystery 
(see 'Alcx anded |2005| for a review; see also IPerets et al.l 
i2006 and references therein.). Outside this cluster, at 1- 
13 arcsec (0.04-0.5 pc), is at least one disk of young stars, 
with a stellar population distinctly younger t han the S- 
stars (|Genzel et al. I I2OO3I : iPaumard et al.|[2006l ). The spec- 
tral identification of > 30 post main sequence blue super- 
giants and Wolf-Rayet stars con strains the age of the disk 
to < 8Myr (|Paumard et aklboOfil) . T he stars have top-heavy 
initial mass function with /3 » 1.35 ([N ayakshin et al."2006|; 
iPaumard et al. I I2OO6I : see also iNavakshi n fc Sunvaov 200^). 
Within the same observations of the disk, IPaumard et al.l 
( 200 6) found the density of OB-type stars to fall off more 
steeply than the observed cusp of stars, with no positive 
detections outside of 0.5 pc. 

Dispersed throughout the stars there should be a group 
of stella r -mass BHs brought there through mass segregation 
(iMorrid Il993l : iMiralda-Escude fc Gould [2O0O). Assuming 
that all the stars in the ste llar cusp a round Sgr A* formed 
10 Gyr ago, ,Miralda-Escude fc Gouldl ((2OOQ ') predicted that 
there should be ~ 2.5 x 10* BHs within the central pc of 
the Milky Way. In their analysis, iMiralda-Escude fc Gouldl 
(2000) estimated that 1.6% of the stellar mass is in BHs 
given an approximately Salpeter MF. All BHs within 5pc 
relax to the central pc through mass segregation, where they 
form an a — 1. 75 — 2.0 cusp. More detailed stu dies support 
this argument. iHopman fc Alexanderl (|2006lJ ) numerically 
solved the time-dependant Boltzmann equations for a four 
mass model near Sgr A* with 10 Af© BHs, 1.4 M© neutron 
stars, 1 Mq stars, and 0.6 Mq white dwarfs. They found that 
there should be f» 1.8 x 10^ BHs and 3 x 10** Mq of stars 
within 0.1 pc of Sgr A* in an obh ~ 2.0 , a* ~ 1.4 cusp re- 
spectively. In addition, large A'^ Monte-Carlo simulations of 
relaxation and stellar evolution in the Galactic centre foun d 
the formation of a similar cusp of BHs (|Freitag et al.|[2006l ). 
However, in these simulations the density of BHs at 0.1 pc 
was found to be ten times higher a nd the density of stars 
t o be a bout five times lower than in lHopman fc AlexandeJ 
(|2006bl ). Very few observations constrain number of stellar- 
mass BHs considered here, however X-ray observations seem 
to place an upper limit of ^ 40, 000 BHs in the inner parsec 
(|Deegan fc Navakshinll2007l ). 
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3 THEORY AND ASSUMPTIONS 

To calculate the rate at which stars are ejected from the 
Galactic centre through t heir scattering o f f BH s there, we 
generalise the analysis of IYu fc Tremain3 (|2003|) to mu lti- 
mass systems (in a form similar to that of lHenoDlll969l ). In 
our calculations we assume that Sgr A* ha s a mass Mambh = 
3.5 X 10^ Mq and is at a distance of 8 kpc (|Ghez et al.ll2005l : 
lEisenhauer et al.ll2005l l. 



3.1 Two— body scattering 

We would like to identify the conditions under which a scat- 
tering between a BH and a star would result in the ejection 
of a HVS. In each scattering, a st ar of mass m and velocity 
V undergoes a change in velocity (|Binnev fc Tremaindl 19871 . 
Eqs. 7- 10a - 10b) 



Sv = 



2m' bw^ 
G{m' + m)2 

2m.'w 



1 + 



G2(m' + m)2 

b^w^ 
G2(m' + m)- 



(1) 



where b is the impact parameter of the encounter, w — \w\ — 
\v ~ v'\ is the relative velocity at infinity, m' is the mass of 
the BH, and v' is the BH's velocity. The direction of Sv 
is determined by w and w±, the unit vectors along and 
perpendicular to w, respectively. 

In some instances, the BH may tidally distort (or even 
disrupt) the incoming star, causing the BH-star system to 
lose energy. At closest approach, the stars will have a rela- 
tive velocity of ^y w'^ + 2G{m + m') /bmin, where &min is the 
distance of closest approach between the star and BH. In 
order to be conservative, we assume that the star is com- 
posed of two halves, each with mass m/2, separated by a 
distance of two stellar radii, 2_R* . This is chosen in order to 
have the energy loss diverge when bmin ~ R*- Then, in the 
impulse approximation, the total amount of energy lost is 
approximately 



AE: 



2G^mm'^Ri 



(^2 + 2G(m -f m')/femin)(6Ln - -R')' 



(2) 



Thus, the final relative velocity between the star and BH is 
reduced to. 



= ^w^ - 2AE/r. 



(3) 



If we now assume that the final relative velocity is in the 
same direction as in Equation ([T}, we get the final change 
in velocity. 



<5'i; — 



2m'bwf'uP 
G{m' + mY 



b'w" 



m' + m 



w — Wf + 2w / I 1 + 



G2(m'-f m)2 

f,2,„4 \ -1 



b'w* 



w. (4) 



Equation @ correctly reduces to Equation ^ as Wf ap- 
proaches w. 

We are seeking circumstances under which the star's 
velocity at infinity is greater than some threshold ejection 
speed Vej, 



2 I , c |2 2Gi\fsmbli ^ 2 
Voc = \v+dV >Vci: 



(5) 



where r is the distance from Sgr A* at which the encounter 
occurred. For the distance of 55 kpc, at which many HVSs 
have been observed, the star would have a velocit y «55 ~ 
[uL - (800kms-i)Y/2 (|Carlberg fc Innanen|[l987l ). 

For our calculations we determine Rt by fitting a bro- 
ken power-law to t he solar metallicity stellar models of 
ISchaller et all (|l992l ). and do not consider ejections where 
the star is tidally disrupted by the BH (i.e. AE > Gm? /R,). 
For a typical star of 1 Mq at a distance of .01 pc, this con- 
dition would give a minimum closest approach of femin ~ 
2.07ii*. At .001 pc, the star and BH could get as close as 
1.57ii* before tidal disruption. 



3.2 Total rate 

The population of stars near Sgr A* can be described by the 
seven-dimensional distribution function (DF) /, (r,D,m), 
where m is the star's mass, v is the star's velocity, and 
r is the star's position relative to Sgr A* which is lo- 
cated at r = 0. Similarly, we describe the DF of BHs by 
fw{r' ,v' ,m'). 

In our analysis we restrict our attention to stars scat- 
tering off BHs only, since in star-star scattering, physical 
collisions between the stars limit the total production rate of 
HVSs (Yu & Tremaine 2003). Thus, the probability for a test 
star to encounter a BH at an ii npact pararaeter b within an 
infinitesimal time interval d£ is (|Henonll960l : lYu fc Tremaind 
I2OO3I ) 

T{r,v)At = At j bdb J wA^v' j d* j Am' fBnir ,v' ,m!), 

(6) 

where m' is the BH mass, and 'I' is the angle between the 
(i), i;')-plane and (v — u', <5w)-plane. The integrated rate at 
which all stars undergo such encounters in a small volume. 



d r, around position r is 



7^(r) 



dm ft{r,v, m) V{r, v). 



(7) 



Finally, the total rate of ejecting stars with velocity ^ Uoj is 

(8) 



'^^"J = / dV7^(r) 



At 

from some inner value rmin (which depends on the distribu- 
tion function of the stars, see, H3.3P to 0.1 pc, where we limit 
the integration over b, w, and 'I' such that Voa > Vcj, as in 
Eq. ([5]), and exclude all encounters that tidally disrupt the 
star {AE > Gm^/R,). 

In order to find the total rate of HVSs we integrate 
Eq. ((8)1 using a multidimensional Monte-Carlo numerical in- 
tegrator. Because the integrand of Eq. |S]is not continuous, 
other integration routines are unsuitable for our calcula- 
tions. We integrate all of our results until they converge 
to a residual statistical error of < 10 %. 



3.3 Distribution functions 

Unfortunately, there are no observations to constrain the ex- 
pected number or distribut ion of BHs near Sgr A* (see also, 
iDeeean fc NavakshinllioOTh . For simplicity, we assume that 
all BHs have a mass of 10 M©, are distributed isotropically, 
and follow an qbh cusp density profile so that 

/bh oc E"^^-^-"5{m' 



10M(7 



(9) 



4 O'Leary & Loeb 



Table 1. Stellar DF Models and Results. The first four columns describe the different stellar distribution functions used in our analysis. 
/3 is the slope of the MF, which is normalised by the maximum mass (column 3) and minimum mass (column 4) of the stars, as well as 
the total number of stars, as described in ^ 13.31 The remaining columns list the total integrated rate of HVSs (vca > 1000 km s~^) for 
each model. From left to right the rates correspond to m' = 10 Mq, m' = 15 Mq both with r^^^ = 0.001 pc, and finally m' = 15 Mq 
with r„jin = 0. Our results are described in more detail in § |4] 



Model 


/3 


Maximum Mass 


Minimum Mass 




HVS Rate 












m' = 10 Mq 


m' = 15 Mq 


^■min = 






(Mq) 


(Mq) 


(yr-i) 


(yr-i) 


(yr-i) 


la 


4.85 


20 


3 


3.5 X 10-s 


1.5 X 10"* 


8.2 X 10"* 


lb 


3.85 


20 


3 


2.0 X 10-3 


7.8 X 10-3 


3.9 X 10"* 


2a 


2.35 


3 


0.5 


3.3 X 10-'^ 


9.8 X lO"'^ 


3.6 X 10-<5 


2b 


2.35, 4.85 


3 


0.5 


4.4 X 10-'^ 


1.3 X 10-6 


4.8 X 10-8 



where E = (GMsmbh/r — ^'^ /2) is the negative specific en- 
ergy of a BH and obh ~ 2. Equation ((9)1 is normalised 
by integrating Equation [9] over \v'\ < y/2GM/r and r so 
that there are Nbus = 1800 BHs within 0.1 pc of Sgr A* 
consis tent with the calculations of iMiralda-Escude fc Gouldl 
(|2000ll and iHopman fc Alexande il l|2006bl ). In addition, we 



look at a system of A^bhs = 1800 more massive BHs, 
m' = 15 Mq, in order to determine the importance of the 
BH MF on the rate of HVSs. Because the density power-law 
assumed here obh = 2 is steeper than the ~ 1.5 slope de- 
rived for a "drain limite d" population of BHs (see Eq. 1 of 
[Alexander fc Livid liooi ) . the innermost regions of the in- 
tegration (interior to < lO"* — .001 pc) will have a number 
density that exceeds this theoretical limit. Therefore, at the 
end of § |4] we also look at an approximately drain limited 
sample of BHs. 

The stellar DF, on the other hand, should be much shal- 
lower than the BH's within ^ 0.1 pc of Sgr A*. Obse r vations 
as well as the calculations of iHopinan fc A lexande j 1 2006bl ) 
suggest a* = 1.4 in this region IjGenzel et al. 20031 ) . This 
gives 



/. oc E 



-0.1 dn 
dm ' 



(10) 



where dn/dm oc m is the current MF of stars. We trun- 
cate Eq. (fTO)l so that /* = when 



E > G{m' + m)/R,. 



(11) 



Stars would need to undergo physical collisions with the 
BHs, which are rare a nd disruptive, in orde r to relax to a 
higher specific energy (jBahcall fc Wolj|l976l ). For a typical 
1 Mq star, this corresponds to an orbit with minimum semi- 
major axis « .007 pc, re markably close to the semimajor axis 
of S2 (SO-2) of .005 pc (|Ghez et al.ll2005l : lEisenhauer et al.1 
l2005l V Because we assume that there are no stars with a 
tighter orbit, interior to the region the number density of 
stars is much shallower, and follows an a* ~ .5 power-law. 
In addition to the above constraint, we explicitly exclude all 
orbits which would tidally disrupt the star. 



'closest 



M„ 



1/3 



(12) 



where rdosost is the star's closest approach to the SMBH. 
Although Equations [11] fc [12] should completely determine 
the innermost radius of the stars, rmin, we also set /* = for 
r < rmin. Except where explicitly mentioned in our results, 
we conservatively set rmin = 0.001 pc. This is the innermost 



resolved radius of the simulations of iFreitag et al. I l|2006l ). 
where the authors found no significant change in the DF 
caused by stellar collisions. 

For the stellar DF, we use a total of four models, varying 
both /3 and the normalisation. We do this both to constrain 
the total number of HVSs as well test how sensitive the 
rate is on the model. For clarity of our results, as well as to 
satisfy different observational constraints, we separate the 
stars into two distinct populations. 

For stars with mass m > 3 Mq , we are constrained 
by the total number of observed S-stars, as well as the 
star formation history in the region. There are 0(10) S- 
stars with mass m > 5 Mq within .01 pc of Sgr A* 
(Ghcz et al. 2005; Eisenhauer et al.l [20051 ). We assume that 
this reflects approximately a steady state and normalise our 
models (la-b) to 10 such stars subject to the above con- 
straints (i.e., Eqs. [TT] fc I12p . To have a complete descrip- 
tion we also need to know the slope of the current MF of 
stars, p. Unfortunately, this is unconstrained for stars with 
m < 5 Mq . However observations suggest that all of the 
stars in the region have formed mostly cont inuously in time 
dAlexander fc Sternberg 19991 : [Genzel et al.| l2003). For com- 
pleteness, we look at two different values of /3. If the stars 
formed with a standard Salpeter MF continuously through- 
out time, we would expect dm/dn oc m~^'^^ x tm, where tm 
is the lifetime of the stars. For m < 10 Mq, tm. oc m~^'"', 
therefore for Model la, we set /3 = 2.35 -I- 2.5 = 4.85. 
In Model lb, we set (3 = 1.35 -1-2.5 = 3.85, consistent 
with a fla tter MF as observed in the young disks of mas- 
sive stars ([Navakshin et al.[[2006l : [Paumard et al.[[2006l ). We 
have conflrmed that the total amount of background light 
contributed by the unresolved stars is very near, but does 
not exceed the best esti mates from the obse rved background 
light in the inner cusp l|Schodel et al.[[2007i ). 

For stars with mass 0.5 Mq < m < 3Mq, we use two 
different models. For Model 2a, we use a Salpeter MF with 
P = 2.35. For Model 2b, we use a broken-power law with 
(3 = 2.35 for 0.5 Mq < m < 1 Mq, and /S = 4.85 for 1 Mq < 
m < 3Mq, roughly consistent with a population of stars 
that formed continuously over the last 10 Gyr. In both cases, 
we normalise Eq. (|10p by requiring that the total mass in 
stars for r < 0.1 pc be equal to 6 x 10* Mq ([Genzel et al.[ 
[2003h . 
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4 RESULTS 

In our results, we consider a star to be a HVS when Voo > 
lOOOkms"^ {vcj = 1000 kms"^). We list our results alone 
with a summary of each Model's MF in Table [1] In Figure [T] 
we show how the total ejection rate depends on the minimum 
ejection velocity, «cj, with r-min = 0.001 pc. For small ejec- 
lOOkms"^) the rate decreases roughly 



tion velocities {vc 
as a power- law o 



independent of the stellar MF or 



mass of the BHs. However, for velocities Vcj > 800kms~ 
tidal dissipation as well as physical collisions between stars 
and BHs begin to suppress the rate of HVSs. In our Model 
la, with mniin = SMq and 10 BHs, the ejection rate 
decreases rapidly from 3.5 x 10~^ yr~^ for Vej = 1000 km s~^ 
(with an observed velocity of V55 > 600kms~^ at 55kpc 
from Sgr A*) down to 6.3 x 10"" yr'^ for Uej = 1500 km s"^ 
iv55 > 1300 kms^^). The velocity distribution of lower mass 
HVSs (~ .5 Mq) shows a slightly weaker change in the 
rate of HVSs in the high velocity range, as can be seen 
in Figure [T] Thus, the observed distribution of HVSs pro- 
duced from this mechanism should be truncated sharply 
around ~ 2000 km s~^, in contrast to tidally disrupted bi- 
naries which easily eject HVSs with velocities exceedin g 
4000 kms"^ (|Ginsburg fc Loebll2006l : iBromlev et al.ll2006l ). 

Figure [2] shows the dependence of the rate of HVSs on 
the stars' masses. The rate of HVSs is plotted in 0.5 Mq 
intervals, so that the total rate of a model is the sum over 
all mass bins. For a cluster of single mass BHs, the mass 
distribution of HVSs produced by this mechanism is steeper 
than the stellar MF by ~ 1 for the low mass stars, and 
~ 2 for higher mass stars. A steeper slope is expected since 
more massive stars will receive a much smaller kick than 
less massive stars (~(x m'/(m + m')). However, it remains 
unclear if these slopes will hold for cluster of BHs with varied 
masses, where a few of the most massive BHs may eject the 
majority of HVSs. The steep slope of the HVS MF makes it 
clear that very massive HVSs (m ~ 10 Mq) can not originate 
from this mechanism, unless there were significantly more 
stars in the past. There is also a clear discontinuity in the 
rate of HVSs at 3 Mq . This just reflects the clear overdensity 
of high mass stars in this region. 

The ejection rate also depends very strongly on the ex- 
tent to which stars populate the innermost regions of the 
cusp, as shown in Figure [3] This is despite the constraints 
on the stellar DF (c.f. Eqs. [11] & [12} , which limit both the 
semimajor axis of stars to ~ 0.01 pc and the maximum ec- 
centricity of the orbit. For rmin ^ 0.001 pc the total rate 
increases with a shallow power-law index versus inner ra- 
dius, doubling its value by 10~^ pc. We caution, however, 
that interior to 0.001 pc , stellar collisions may flatten the 
cusp of low mass stars l|Murphv et al.l Il99ll : iFreitag et al.l 
I2OO6I ) and cause the rate of HVSs to increase more slowly 
than shown in Figure[3l For the high mass stars in the region 
(Models la-b), the stars are less likely to have physical col- 
lisions in their lifetimes, and may populate the entire phase 
space aUowed by Eqs. ([Til ^ pS)) . 

Smce IBrown et all (|2006al ) initiated a targeted search 
for HVSs in the Galactic halo, it is useful to forecast the 
number of observable HVSs that would originate from en- 
counters with stellar-mass BHs. Current observations indi- 
cate that there are ~ 33 ± 17 hypervelocity B-stars with 
masses between 3-5 Mq (W. Brown, private communica- 
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Figure 1. The total rate of stars with ejected velocities v^c > foj 
versus v^j for r^i^ = 0.001 pc. The coloured curves from top to 
bottom are for Models 2a (green), 2b (magenta), la (black) and 
lb (red). The dashed curves are for the same stellar MFs but for 
a system of 15 Mq BHs instead of 10 Mq. The top axis is labelled 
to show the corresponding velocity at a galactocentric radius of 
55kpc, D55. 
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Figure 2. The total rate of HVSs (v^j = 1000kms"\) in 0.5Mq 
intervals for Tmin = 0.001 pc. The colours correspond to the same 
models as in Figure[T] The top and bottom black 'dotted' lines are 
for a power-law with slope -3, and -7 respectively for comparison. 
The discontinuity of the HVS rate at 3 Mq is a result of our two 
different normalisations for Models 1 and 2. 
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tion) . Can scattered stars alone explain the abundance of the 
observed HVSs? The answer is positive but not under the 
most conservative conditions assumed thus far, which pre- 
dict only ~ l(iVBH/1800) such HVSs, for r^in = 0.001 pc. 
If the masses of the BHs in the cluster were more typically 
~ 15 M0, as might be expecte d from mass seRregatio n as 
well as an old population of BHs (|Belczvnski et al.l2004 ) , the 
number HVSs in our Galaxy would be four times larger. This 
is still lower than the number of observed HVSs. However, if 
we assume the most massive stars populate the entire phase 
space within the constraints of Eqs. (|11[) & (|12|l , we expect 
there to be ~ 30 HVSs in our Galaxy. If we relax our defini- 
tion of a HVS and set Ve\ = 900kms~^, the ejection velo city 
of HVS7 (|Carlberg fc Innanenll987l : lBrown et aLll2006bl ). we 
get ~ 50 HVSs. The assumption that the stars fill the entire 
phase space is in best agreement with the observ ations of the 
S-sta rs, which appear to be a relaxed system (iGhez et al.l 
I2OO5I ). By integrating the DF in this region, subject to the 
above constraints, we find that there should be < 0.3 stars 
interior to 0.001 pc. In addition. Equation limits the 
semimajor axis of any star to ~ O.Olpc, close to the semi- 
major axis of S2 (SO-2), ~ 0.005 pc. Thus, the stars that are 
ejected from this region are ones on eccentric orbits which 
bring them into this region, similar t o S14 (SO-16), wh ich has 
a closest approach of 0.0002 pc l|Ghez et ahlbOOSl ). This 
is also cons istent with the re sonant relaxation timescales in 
the re gion (jRauch fc Tremaine 1996; Hopman fc Alexander 
l2006al . however, note our caution in §[5}. We also note, how- 
ever, that we are still being conservative in our calculations. 
We have excluded all collisional encounters, which at the 
highest relative velocities may not actually disrupt the star 
or have significant tidal energy dissipation. If stars and BHs 
do populate the inner 10~^ pc, any collisions would be at ve- 
locities far greater than the escape velocity of the star. Our 
calculations suggest that such encounters may contribute a 
factor of a few more HVSs, and may offset any depletion of 
stars and BHs in the region. Since any such collisions would 
be brief and not very luminous (unless coalescence follows at 
low impact speeds) , it would be difficult to identify star-BH 
collision events in external Galactic nuclei. 

Most recently, iBrown et al.l l|2007h have detected a sig- 
nificant number of marginally bound HVSs in the galactic 
halo. They found that these stars outnumber the unbound 
population by a factor of ~ 2. Such a ratio of bound to un- 
bound HVSs is consistent with the scenario presented here, 
as well as the tidally ejected binary scenario iHillsl jT988) to 
within the statistical uncertainties (iBrown et al. [20071 ) . If 
we look at the ratio of the rates of stars ejected with veloc- 
ities between 845kms~^ (1)55 = 275kms~^) and 920kms~^ 
('i'ss = 450kms~^), and stars ejected with higher velocities, 
we find that there should be ~ 2.6 times more bound HVSs 
than unbound HVSs in a volume limited survey. However, 
this analysis assumes that all such HVSs would be observ- 
able, and does not account f or any selection effects of the 
survey bv lBrown et~aLl (|2007h . With more HVS discoveries, 
a more detailed statistical analysis may be able to determine 
the origin of the HVSs. 

As we pointed out in § 13.31 the density profile 
of BHs assumed thus f ar will exceed the "drain limit" 
l| Alexander fc Livid l200^ ) in the innermost regions of inte- 
gration. We have therefore approximated the drain limit for 
both 10 M0 and ISM© BHs by using an qbh = 1.7 power- 
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Figure 3. The HVS ejection rate {v^j = 1000 km s"^) versus the 
minimum radius of integration r^i^. The lines are coloured as in 
Fig.[T] Interior to .001 pc, stellar col lisions may flatten the cusp 
of low mass stars llMurphv et al.11199 1: Frcitag ct al. 2006) and 
cause the rate of HVSs to increase more slowly than shown in 
Figure|3] For the high mass stars in the region (Models la-b), the 
stars are less likely to have physical collisions in their lifetimes, 
and may populate the entire phase space allowed by Eqs. mil l &c 

m- 

law, whic h we found to b e the best fit to the Equation 1 of 
Alexand er fc Livid (|2004l ) between 0.001 and 0.1 pc. We as- 
sumed that, inside of 0.1 pc, there are A'bh ~ 12000 10 M© 
BHs and « 5500 15 Mq BHs. We have found that this in- 
creases the total rates calculated above by a factor of ~ 4 for 
the 10 Mq BHs and - 2 for the 15 Mq BHs. This is true, 
even if the BHs are drain limited only out to ~ 0.01 pc. 
Thus, for a drain limited population of 10 Mq BHs (inte- 
rior to O.Olpc), there would be ~ 40 hypervelocity B-stars 
produced from BH-star scattering, with rmin = 0. 



5 SUMMARY AND DISCUSSION 

We constructed a variety of models for the distribution of 
stars and BHs in the innermost 0.1 pc of the Galactic cen- 
tre. Assuming that the cusp and disk of stars we observe 
today represent a steady state over the past ~ 100 Myr, we 
calculated the rate at which stars will scatter off BHs and 
populate the Milky Way halo. We showed that the total 
ejection rate of HVSs can be comparable to that produced 
by the tidal disruption of binaries, and may exceed it for 
intermediate mass stars. Our results are consistent with the 
total number of ejected stars from Fokker-Plank simulations 
by .Freitag et aL (200Q). Assuming that the B-st ars are fully 
relax ed as is observed with the younger S-Stars (|Ghez et al.l 
l2005h . we can account for mos t, if not all, of the stars ob- 
served bv lBrown et al.l l|2006ah . We demonstrated that the 
velocity distribution as well as the mass distribution of HVSs 
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should be truncated at high values for the BH scattering pro- 
cess compared to binary disruption events (see Figs. 1 and 
2). Better statistics of HVS detections could therefore deter- 
mine the relative significance of these two plausible channels. 

In our study, the ejected stars originate from the in- 
ner ~ 0.1 pc near Sgr A*, whereas the tidally disrupted 
binaries that produce HVSs originally come from > 2 pc 
||Yu fc Tremainell2003l : iPerets et al 20061. The o bservations 
of the young disk of stars ( Paumar d ct al.| |2006l ). as well as 
the cluster of S-stars, suggest that the population of stars 
at 0.1 pc from Sgr A* should be younger and more massive 
than at 2 pc. This is furt her supported by the steep drop in 
B-stars outside of 0.5 pc (jPaumard et al.ll2006l '). In addition, 
binaries with B-stars (tms ^ 300 Myr) may also not be fully 
relaxed as has been assumed in many previous calculations, 
and therefore the diffusion rate of the most massive binaries 
into the loss-cone may be lower by orders of magnitudes. 

In the survey of iBrown et~aLl l|2006bl ), the observed 
HVSs have the same colour as blue horizontal branch stars 
(BHBs) and have yet to be distinguished from B-stars with 
high resolution spectroscopy. The calculations shown here 
suggest that low mass stars, such as BHBs and their pro- 
genitors, can become HVSs with high efficiency if they are in 
the inner 0.1 pc of Sgr A*. Although tidal disruption rates 
may be enhanced by massive perturbers, BHBs must un- 
dergo significant mass loss, causing them to widen or even 
disrupt, and therefore are un likely to be found in the tight 
binaries which are disrupted (jPerets et al.ll200^ . The spec- 
troscopic identification of one hypervelocity BHB star would 
be strong evidence in support of the mechanism presented 
here, or an inspiralling massi ve BH l|Yu fc TremaindbOOj : 
lLevin|[2005l : iBaumgardt et al.ll2006l ). However, the expected 
rate of hypervelocity BHBs depends on uncertain details of 
stellar evolution and goes beyond the scope of this work. A 
population of BHBs ejected through this mechanism may 
be spun up through repea ted encounters in the BH cusp 
l| Alexander fc Kumaiil200ll ). as well as during the final en- 
counter which ejects them (H. Perets, private communica- 
tion). Therefore, a rotating BHB star may appear to be more 
like a B-star without det ailed spectroscopic evid ence. In fact, 
the HVS HE 0437-5439 (|Edelmann et all 120051 ') . should not 
be ruled out as being a BHB star based on its required spin 
value alone. 

There is a curious link between the number of HVSs 
observed by iBrown et al.1 l|2006bD and the S-stars ob- 
served orbiting Sgr A*. Some S-stars may be the for- 
mer companions to the HVSs from tidally dis r upted 
binaries (Gould fc QuiUenI l2003l : iGinsburg fc Loebl l2006l : 
iPerets e~al.ll2006i ). Interestingly, there exists a similar con- 
nection between HVSs scattered off of the BHs and the S- 
stars as well, although perhaps not one-to-one. For every 
strong encounter that produces a HVS, there is likely an- 
other encounter which can bring the star closer to Sgr A*, 
perhaps kic king out the BH instead, s imilar to the scenario 
proposed bv [Alexander fc Livid (|2004h . In this scenario, the 
scattering of a small fraction of young stars from many pre- 
vious disks into the inner 0.04 pc may result in a population 
similar to the S-stars (O'Leary, R. M. fc Loeb, A. 2007, in 
preparation). However, we should note that if the S-star are 
in fact the remnants of tidally disrupted binaries, then this 
mechanism will only be secondary, as it is only ejects a small 
fraction of such stars. 



There is obviously considerable uncertainty in our mod- 
els. The rate of diffusion of stars both close to and far from 
Sgr A* into eccentric orbits is important to understanding 
both the mechanisms as well as the source of the observed 
B-type HVSs in the Galactic halo. In both cases, it is most 
likely that the stars formed on relatively circular orbits, 
whether in a disk around Sgr A* or in an inspiralling clus- 
ter. Lar ge scale simulations similar to those already done 
bv iFreitag et al.l (|2006| ) can help resolve the uncertainties of 
relaxation, and with modifications, may be able to account 
for the conditions of continuous star formation over long pe- 
riods of time. For low mass stars, physical collisions may 
deplete the total number of stars, and reduce the rate. In 
our discuss ion we have not considered the effects of resonant 
relaxa tion (|Rauch fc Tremainelll996l : iHopman fc Alexander] 
I -'I II H.. J ) near Sgr A*, which may have two counteracting ef- 
fects on our rate calculation. Resonant relaxation may flat- 
ten the cusp and deplete the number density of stars and 
BHs in the innermost 0.01 pc; at the same time, it may also 
drive more massive stars into the same region producing 
more B-type HVSs (|Hopman fc AlexandeIj|2006a^ . 

In our analysis, we also neglected the migration 
of massive objects n ear S gr A*. If, as suggested by 
IPortegies Zwart et al.l (|2006l ). many intermediate mass BHs 
(IMBHs with masses > 10^ M©) populate the inner pc of 
the Galactic centre and merge with Sgr A* every ~ 10^- 
10* yr, then the cusp of stellar-mass BHs and stars would 
not regenerate f ast enough to produce HVSs through BH- 
star encounters (|Baumgardt et al]|2006l ) but i nstead coul d 
produce them through IMBH-star encounters (|Levinll2005l '). 
BHs with masses between most stellar mass BHs and IMBHs 
(e.g. ~ 50 Mp) ) , may also form through standard binary evo- 
lution (|Belczvnski et af]|2004h . Such BHs could considerably 
increase the rate of HVSs, even if they are relatively rare. 
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